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ABSTRACT 

The recent long-duration GRB 100316D, associated with supernova SN 2010bh and detected by 
Swift, is one of the nearest GRB-SNe ever observed (z = 0.059). This provides us with a unique op- 
portunity to study the explosion environment on ~ kpc scale in relation to the host galaxy complex. 
Here we present spatially-resolved spectrophotometry of the host galaxy, focusing on both the explo- 
sion site and the brightest star-forming regions. Using these data, we extract the spatial profiles of the 
relevant emission features (Ha, H/3, [OIIIJA5007, and [NII]A6584), and use these profiles to examine 
variations in metallicity and star formation rate as a function of position in the host galaxy. We con- 
clude that GRB 100316D/SN2010bh occurred in a low-metallicity host galaxy, and that the GRB-SN 
explosion site corresponds to the region with the lowest metallicity and highest star formation rate 
sampled by our observations. 
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1. INTRODUCTION 

Long-duration gamma-ray bursts (LGRBs) arise from 
the core-collapse deaths of some massive stars. The pri- 
mary line of evidence for this association comes from 
the detection of supernova signatures following the GRB. 
These include photometric re-brightening with the mag- 
nitude, color, shape, and timescale expected for a super- 
nova (e.g. Bloom ct al. 1999, Zeh et al. 2004, Soderberg 
et al. 2006, Woosley & Bloom 2006, Cobb et al. 2010), 
and more importantly from the spectroscopic detection of 
broad-lined Type Ic supernovae (SNe) in several LGRBs 
at z < 0.3 (e.g., Galama et al. 1998, Stanek et al. 2003, 
Malesani et al. 2004, Modjaz et al. 2006). The connec- 
tion with Type Ic SNe is extremely important, since these 
events are thought to be caused by the collapse of mas- 
sive stars that have shed their outer hydrogen and helium 
envelopes (see Filippenko 1997). Similarly, the broad ab- 
sorption features indicate the presence of large ejecta ve- 
locities, and perhaps energies (e.g. Galama et al. 1998, 
Pian et al. 2006, Modjaz et al. 2008). 

The most recent addition to this class of spectroscop- 
ically identified GRB-SNe is GRB 100316D/SN 2010bh 
at z = 0.059 (Stamatikos et al. 2010, Vergani et al. 2010, 
Bufano et al. 2011). In Chornock et al. (2011) we re- 
ported the spectroscopic discovery of the accompanying 
SN 2010bh, classified it as a broad-lined Type Ic SN 
based on optical and near-infrared spectra, and delin- 
eated the evolution of its photospheric velocity. We also 
demonstrated that the host galaxy has a low luminosity, 
L w 0.1 L», and a low metallicity, Z < OAZq, similar to 
other GRB-SNe (e.g. Levesque et al. 2010a). In addition, 
Starling et al. (2011) presented Hubble Space Telescope 
observations of the host galaxy, which revealed a complex 
morphology, possibly suggestive of recent merger activ- 
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ity 

The metallicities of the GRB 100316D/SN 2010bh host 
galaxy and explosion site are of particular interest. In 
recent years, a connection has been suggested between 
LGRBs and low-metallicity host environments. LGRBs, 
associated with star-forming host environments, are ob- 
served in fainter and more irregular galaxies than their 
core-collapse SN counterparts (e.g. Fruchter et al. 2006, 
Wainwright et al. 2007). Their host galaxies, on av- 
erage, fall below the luminosity-mctallicity and mass- 
metallicity relations for star forming galaxies out to z ~ 1 
(e.g. Modjaz et al. 2008, Kocevski et al. 2009, Levesque 
et al. 2010a, b). However, the physical mechanism driv- 
ing this apparent metallicity trend is still poorly un- 
derstood. LGRB host environments do not appear to 
adhere to any low-metallicity cut-off, with several high- 
metallicity host galaxies and explosion sites noted in the 
literature (Graham et al. 2009, Levesque et al. 2010b, c). 
Similarly, the relativistic supernova SN 2009bb, found 
to have a comparable ejecta velocity to GRB-SNe but 
with no accompanying LGRB, was also observed in a 
high-metallicity host environment (Soderberg et al. 2010, 
Levesque et al. 2010d). There is also no apparent cor- 
relation between host metallicity and gamma-ray energy 
release for LGRBs (Levesque et al. 2010e), a result that is 
at odds with previous predictions of the progenitor mod- 
els (MacFadyen & Woosley 1999). The low-metallicity 
trend could also be an artifact of these events young pro- 
genitor ages, the star formation histories of their hosts, 
or selection effects related to their star formation rates 
(e.g. Bloom et al. 2002, Berger et al. 2007, Kocevski & 
West 2011). 

One limitation of the aforementioned studies is their re- 
liance on global metallicities. For the majority of LGRB 
hosts at z > 0.3, pinpointing the LGRB explosion site 
and acquiring site-specific spectra within the small, faint 
host galaxies is a difficult proposition. The nearby GRB- 
SNe therefore offer a unique opportunity to probe the 
metallicities of their host environments in detail. For 
these spatially-resolved hosts we can determine metallic- 
ities at the LGRB host site as well as the surrounding 
star-forming regions of the galaxy, allowing us to pin- 
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point the precise environments that produce GRB-SNe 
and place these host sites in context with their global 
host galaxy properties. 

Here we present spatially-resolved spectroscopy of the 
GRB 100316D/SN 2010bh host galaxy. These observa- 
tions allow us to place the properties of the explosion site 
in the context of the broad host properties. We discuss 
the observations and describe the analysis technique in 
§2. We derive spatially-resolved line fluxes, line ratios, 
metallicities, and star-formation rates (SFRs) across the 
host galaxy, with a particular focus on the GRB-SN ex- 
plosion site in §3. Finally, we consider the implications 
that these results for our understanding of LGRB host 
galaxies and the progenitor population in §4. Through- 
out the paper we use the standard cosmological parame- 
ters H = 71 km s" 1 Mpc -1 , TO = 0.27, and fl A = 0.73. 

2. OBSERVATIONS AND REDUCTION 

On 2010 May 8 UT we obtained a 2400 s spectrum of 
the GRB100316D/SN2010bh explosion site (hereafter, 
"site"), using the Low Dispersion Survey Spectrograph 
(LDSS3) on the Magellan/ Clay 6.5-m telescope at Las 
Campanas Observatory, using the 1" center slit; seeing 
was 0.88" (as measured from acquisition images in r) 
and the spectral resolution was 8.5A. A second spectrum 
was obtained on 2010 May 10 UT for a total exposure 
time of 1800 s, with the 1" blue slit aligned along the 
extended bright region of the host galaxy complex (here- 
after, "host"); seeing was 0.62" (as measured from acqui- 
sition images in i) and the spectral resolution was 9.0A. 
The slit positions for the "site" and "host" observations 
are shown in Figure 1 (top left and top right, respec- 
tively). As the airmass was low (~ 1.4— 1.5) the impact 
on the spectra from not observing at the parallactic an- 
gle is expected to be minimal. We also acquired obser- 
vations of the spectrophotometric standard stars EG 131 
and LTT 3864 (Bessell et al. 1999, Hamuy et al. 1994) 
for flux calibration of the "site" and "host" spectra re- 
spectively. All spectra were taken using the using the 
VP H- All grism, and cover the wavelength range 3800- 
7200A. 

To locate the explosion site of GRB 100316D/SN 
2010bh along the slit in the "site" spectrum we astro- 
metrically aligned the r-band acquisition image with an 
r-band image of the GRB-SN from 2010 March 19.99 UT 
obtained with the same instrument. Using 65 objects in 
common to the two images, we find an astrometric accu- 
racy of la — 16 mas in each coordinate. This is signifi- 
cantly better than a single pixel (0.19"), but we expect 
an overall systematic uncertainty of about ±1 pixel due 
to the overall slit positioning. 

We reduced the data using standard routines in IRAI0, 
including bias correction, flatfielding, and cosmic ray re- 
moval. To isolate the host galaxy emission lines, and 
to remove contamination from a foreground star in the 
bright southwestern region of the host complex in our 
"site" observations (see Figure 1), we fit the continuum of 
both objects with a legendre polynomial using the IRAF 
task background and used this fit to subtract the contin- 
uum from our two-dimensional spectra. The results are 
shown in Figure 2. 

2 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under cooperative agreement with the NSF. 



We extracted line profiles along the spatial direc- 
tion for four emission lines in each of our observed 
two-dimensional spectra: H/3, [OIII]A5007, Ha, and 
[NII]A6584. The bright Ha and [OIII]A5007 features were 
traced using the IRAF task apall with an optimal ex- 
traction algorithm that identified and rejected deviant 
pixels. For the weaker [NII]A6584 and H/3 features, the 
line profiles were extracted by applying the robust spa- 
tial profile of the Ha and [OIII]A5007 features, respec- 
tively. Flux calibration was then applied using a sensi- 
tivity function of the standard stars EG 131 and LTT 
3864 with the IRAF tasks setairmass, standard, and 
sensf unc and determining the response of the pixels cor- 
responding to the emission line wavelengths. Finally, the 
flux-calibrated line profiles were corrected for a Galactic 
foreground extinction of E(B — V) = 0.116 (Schlegel et 
al. 1998, Chornock et al. 2010). Using the ratio of the 
Ha and H/3 features, we determined that host extinction 
was negligible across both of our host profiles and re- 
quired no additional correction of the line profiles; this is 
in agreement with the results of Chornock et al. (2010), 
who also measured a negligible level of extinction in the 
explosion site HII region. Our final line profiles are plot- 
ted in Figure 1. 

3. SPATIALLY-RESOLVED ISM PROPERTIES 

We used our extracted profiles of the H/3, [OIII] A5007, 
Ha, and [Nil] A6584 emission lines to construct emis- 
sion line diagnostic ratio profiles for both the "site" and 
"host" observations. Our two diagnostic ratios of inter- 
est are [Nil] A6584/Ha and [OIII] A5007/H/3. The [Nil] 
A6584/Ha ratio is strongly correlated with metallicity 
(Veilleux & Osterbrock 1987, Kewley et al. 2001), due 
to the dependence of the [Nil] flux on primary and sec- 
ondary nitrogen production (e.g. Chiappini et al. 2005, 
Mallery et al. 2007). However, Kewley & Dopita (2002) 
also note that, due to the low ionization potential of [Nil] 
A6584, this ratio is also somewhat sensitive to the ion- 
ization parameter of the host environment Similarly, 
[OIII] A5007/H/3 is sensitive to the hardness of the ion- 
izing radiation field, making it a useful tracer of the ion- 
ization parameter (Baldwin et al. 1981). While this ratio 
can also be double- valued with metallicity, it is far more 
sensitive to the ionization parameter at lower metallic- 
ities (Kewley et al. 2004). Since the ionization param- 
eter itself is inversely dependent on metallicity (Dopita 
et al. 2006), these two ratios ultimately allow us to iso- 
late metallicity as the fundamental physical parameter 
driving the evolution of our spatial profile. Another ad- 
vantage to both of these diagnostic ratios is their insen- 
sitivity to reddening corrections, due to the close prox- 
imity of the emission lines being compared in each case 
(we found that effects from both the host and the fore- 
ground Galactic extinction were negligible in our anal- 
yses). With profiles for both the [Nil] A6584/Ha and 
[OIII] A5007/H/3 diagnostic ratios, we used the Pettini 
& Pagel (2004; hereafter PP04) calibration of the (93iV2 
metallicity diagnostic to construct metallicity profiles for 
the regions of the GRB 100316D host complex covered by 
our "site" and "host" observations. Finally, we use the 

3 Here we define ionization parameter as the maximum velocity 
possible for an ionization front being driven by the local radiation 
field. 
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Ha fluxes in our data to calculate SFR profiles for our 
observations based on the relation of Kennicutt (1998). 

The results of these analyses are plotted in Figure 1 
(bottom). For our calculation of the emission line diag- 
nostic ratios, PP04 metallicity, and SFR, we restricted 
our analyses to regions where all emission features had 
a > 5<T detection; this was generally limited by the flux 
of the [Nil] A6584 line. The data were also smoothed 
with a boxcar average of 5 pixels in the spatial direc- 
tion, corresponding to the seeing conditions during the 
observations. It can be seen from these observations that 
GRB 100316D occurred in a region of the host complex 
with low metallicity and high SFR. From our spatial pro- 
files and the position of the GRB-SN, we find a PP04 
metallicity of log(0/H) + 12 = 8.2 ± 0.1 at the explo- 
sion site. This in excellent agreement with the metal- 
licity determined in Chornock et al. (2011) based on an 
LDSS3 spectrum taken 3.3 days after the explosion at 
the GRB-SN explosion site. From examining both the 
"site" and "host" metallicity profiles, the host galaxy 
has a relatively flat metallicity profile, with a variation 
of approximately 0.3 dex and an average PP04 metal- 
licity of log(0/H) + 12 w 8.3 ± 0.1; these results also 
agree with the metallicity determined by Starling et al. 
(2011) in another region of the host complex. For both 
profiles, the minimum PP04 metallicity corresponds to 
the bright region of the host galaxy associated with the 
GRB 100316D/SN 2010bh explosion site. We also find a 
SFR of ~1.7 M yr" 1 at the host site, the highest SFR in 
the host complex. This is consistent with SFRs for other 
LGRB host galaxies and explosion sites (Levesque et al. 
2010a, 2010b), although it should be noted that this is 
measured over a small ~ 1" x 1" region (corresponding 
to <~1.3 kpc 2 ), whereas in other host studies the SFR is 
integrated over entire galaxies, with a much larger area 
(-10-30 kpc 2 ). 

It should be noted that different strong-line metallicity 
diagnostics can produce drastically different mctallicities 
from the same spectra and emission line fluxes, with dif- 
ferences ranging up to ~0.8 dex (Kewley & Ellison 2008). 
As a result, we restrict our analyses and comparisons to 
the PP04 metallicity diagnostic. For comparisons with 
mctallicities determined from other diagnostics, the con- 
versions of Kewley & Ellison (2008) should be used. 

4. DISCUSSION 

We have examined the host complex of GRB 
100316D/SN 2010bh, constructing spatially-resolved 
emission line and metallicity profiles from spectrophoto- 
metric observations centered on both the explosion site 
and the brightest region of the host. Based on our re- 
sults, we find little variation in metallicity across the host 
galaxy; however, we note that the GRB-SN explosion 
site is localized near a bright region of the host that 
represents a metallicity minimum and SFR maximum 
in our observations. LGRBs have long been associated 
with star-forming regions due to the young ages of their 
progenitors (<10Myr; Woosley et al. 2002). However, 
the lack of substantial variations in the host metallicity 
and the placement of the explosion site in the lowest- 
metallicity region of the host both have important impli- 
cations for future studies of LGRB host galaxies. 

In Figure 3 we compare host galaxy metallicities and 
explosion site metallicities for a small number of GRB 



host galaxies with spatially-resolved abundances. Using 
integral field unit spectroscopy of the dwarf SBc host 
galaxy of GRB 980425/SN 1998bw, Christensen et al. 
(2008) determined PP04 mctallicities at 23 different sites 
across the host galaxy. Similar to the GRB 100316D/SN 
2010bh host, they find a low average metallicity and no 
evidence of any statistically significant abundance gra- 
dient. They also note that the SN explosion site has a 
lower-than-average metallicity (indeed, only the neigh- 
boring Wolf-Rayet-rich region has a lower metallicity in 
their sample), although their measured abundance agrees 
with the overall host metallicity to within the errors. 
In another recent study, Thone et al. (2008) examined 
spatially-resolved ISM properties across the z — 0.0889 
Sbc spiral host galaxy of GRB 060505, finding that the 
GRB site had a significantly lower metallicity than the 
other regions of the host. However, it should be noted 
that the classification of this particular GRB is a matter 
of hot debate; with a relatively short duration (Tqq = 4 
s) and no accompanying supernova detection despite its 
proximity, the nature of the GRB 060505 and its relation 
to other LGRBs remains unclear (Fynbo et al. 2006). Fi- 
nally, Levesque et al. (2010c) determined PP04 metallic- 
ities for both the nucleus and explosion site within the 
GRB 020819 host galaxy. From Figure 3, we can see 
that for all four of these host galaxies, the explosion 
site metallicities show a good general agreement with 
the galaxies' average (or nucleus, in the case of GRB 
020819) metallicities, displaying a small but uniform off- 
set of —0.1 dex. These results suggest that global metal- 
licities can indeed be treated as representative proxies 
for LGRB host environment metallicities in more dis- 
tant hosts, where spatially-resolved studies are not al- 
ways possible. Interestingly, this also sets the nearby 
LGRB host population apart from the hosts of Type Ibc 
supcrnovae, where large metallicity variations have been 
seen between the explosion sites and host nuclei (Modjaz 
et al. 2011). In addition, we find that for both the GRB 
020819 and GRB 100316D/SN 2010bh host galaxies, our 
highest measured SFR is at the explosion site. Simi- 
larly, Christensen et al. (2008) find a higher luminosity- 
weighted SFR at the GRB 980425/SN 1998bw explosion 
site when compared to the galaxy as a whole or the mean 
luminosity-weighted SFR across the host (though other 
individuals sites do show higher SFRs), and Thone et al. 
(2008) find the highest luminosity-weighted SFR in the 
GRB 060505 host galaxy at the explosion site. 

It is important to recognize that both GRB 980425/SN 
1999bw and GRB 100316D/SN 2010bh represent an un- 
usual class of LGRBs; these events are "sub-luminous" , 
with E %iso values that are much lower that the general 
LGRB population (such events are not even detectable 
at z > 0.3), and show evidence of burst geometries that 
are quasi-spherical, in contrast to their tightly-beamed 
and more luminous counterparts (e.g. Soderbcrg et al. 
2004, Soderberg 2006, Levesque et al. 2010e). This sug- 
gests that these events could be phenomenologically dis- 
tinct from the general LGRB population. However, con- 
tinued studies of these sub-luminous GRB-SNe and the 
ISM properties of their host galaxies could offer impor- 
tant insights into the nature of these events' progenitors 
and their connection with the larger sample of LGRBs. 

Extending these spatially-resolved ISM studies to a 
wider range of nearby LGRB host galaxies would be 
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extremely valuable. The host galaxies of GRB 020903 
(z = 0.251, Ricker et al. 2002), GRB 030329/SN 2003dh 
(z = 0.168, Greiner et al. 2003), and GRB 060218/SN 
2006aj (z = 0.034, Mirabal & Halpern 2006) are all ex- 
cellent candidates for these studies (see also Figure 8 of 
Starling et al. 2011). In addition, while GRB 020903 
and GRB 060218 appear to belong to the same class of 
subluminous GRB/SNe as GRB 100316D/SN 2010bh (a 
supernova association for GRB 020903 was reported in 
Soderberg et al. 2005), GRB 030329/SN 2003dh has a 
notably higher luminosity and is often considered to be 
representative of the higher- redshift "cosmological" class 
of LGRBs (e.g., Stanek et al. 2006); a spatially-resolved 
spectroscopic ISM study of this host galaxy could there- 
fore prove particularly illuminating. Detailed maps of 
metallicity variations in all of these hosts would provide 
important data on the explosion sites of LGRBs, and fur- 
ther clarify how the metallicities of these sites compare 
to overall host metallicities. Additional properties such 
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as ionization parameter, star formation rate, and young 
stellar population age would be useful for characterizing 
the typical explosion environments of LGRBs. Finally, 
stellar population synthesis studies at these sites and 
across the host galaxies would provide a unique oppor- 
tunity to compare the different stellar populations that 
dominate the host light and give rise to the progenitors 
of these rare and intriguing events. 
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Fig. 1. — Top: Acquisition images showing position of the 1" slit for our observations of the GRB 100316D/SN 2010bh explosion site 
(left) and host galaxy (right). The images were taken in the g-band, with contours indicating g-band brightness. Bottom: Emission line 
fluxes, ratios, and metallicity profiles for our spectra of the explosion site (left) and host galaxy (right). From top to bottom, these show 
(a) the Ha and [Nil] A6584 line profiles, (b) the H/3 and [OIII] A5007 line profiles, (c) the [Nil] A6584/Ha ratio, (d) the [OIII] A5007/H3 
ratio, (e) metallicity determined based on the OSN2 diagnostic calibration of Pettini & Pagel (2004), and (f) SFR determined based on the 
Ha diagnostic of Kcnnicutt (1998). On the left, the GRB-SN position is indicated by the solid black vertical line, with errors illustrated 
by the gray region. On the right, the observed metallicity (log(0/H) + f2 ~ 8.2) and SFR (~1.6 M yr _1 ) of the GRB-SN is given by 
the dashed horizontal line in panels (e) and (f), respectively. In panels (a) and (b) flux is in units of 10 — 13 ergs cm -2 s — 1 . 
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Fig. 2. — Background-subtracted 2D spectra of our "site" (left) and "host" (right) spectra, showing the [Nil] A6584 and Ha features (top) 
and the [OIII] A5007.4959 and H/3 features (bottom). In the "site" spectrum, the explosion site of GRB 100316D/SN 2010bh is indicated 
by the dashed black line. The negative residuals in the "site" 2D spectrum are due to subtraction of the bright foreground star (seen in 
Figure 1), and have no impact on the derived line intensities. 
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Fig. 3. — Explosion site metallicities vs. average host metallicities for a sample of four nearby GRB host galaxies. Metallicities for GRB 
980425 and GRB 060505 are taken from Christensen et al. (2008) and Thonc ct al. (2008), respectively, and were determined using the 
Pcttini & Pagel (2004) 03N2 diagnostic calibration. The "host" metallicities here are taken as averages of metallicity measurements made 
throughout the host galaxies. Metallicities for GRB 020819 was taken from Levesque et al. (2010c) and determined using the Pettini & 
Pagel (2004) N2 diagnostic calibration that is better suited for higher-metallicity environments. The "host" metallicity in this case is the 
metallicity of the host galaxy nucleus. Metallicities for GRB 100316D are taken from this work. All four host environments fall ~0.1 
dexbelow the theoretical relation where explosion site metallicity and host metallicity are identical, plotted here as a gray dashed line, 
though this is within the uncertainty of the Pettini & Pagel (2004) diagnostics. 



